arXiv:1502.05597vl [cs.IT] 19 Feb 2015 


On Detection Issues in the SC-based Uplink of a 
MU-MIMO System with a Large Number of BS 

Antennas 


(^)Paulo Torres, (^^Luis Charrua, Antonio Gusmao 
(1,2,3) 1ST - Instituto Superior Tecnico, Universidade de Lisboa, Portugal 
(^)EST - Escola Superior de Tecnologia de Castelo Branco, Portugal 
(^)paulo.torres@ipcb.pt, (^Huischarrua®enautica.pt, (^)gus@ist.utLpt 


Abstract — This paper deals with Single Carrier 
(SC)/Frequency Domain Equalization (FDE) within a Multi- 
User (MU)-Multi-Input Multi-Output (MIMO) system 
where a large number of Base Station (BS) antennas is 
adopted. In this context, either linear or reduced-complexity 
iterative Decision-Feedback (DE) detection techniques 
are considered. Regarding performance evaluation by 
simulation, appropriate semi-analytical methods are 
proposed. 

This paper includes a detailed evaluation of Bit Error 
Rate (BER) performances for uncoded 4-Quadrature Am¬ 
plitude Modulation (4-QAM) schemes and a MU-MIMO 
channel with uncorrelated Rayleigh fading. The accuracy 
of performance results obtained through the semi-analytical 
simulation methods is assessed by means of parallel con¬ 
ventional Monte Carlo simulations, under the assumptions 
of perfect power control and perfect channel estimation. 
The performance results are discussed in detail, with the 
help of selected performance bounds. We emphasize that 
a moderately large number of BS antennas is enough to 
closely approximate the Single-Input Multi-Output (SIMO) 
Matched-Filter Bound (MFB) performance, especially when 
using the suggested low-complexity iterative DE technique, 
which does not require matrix inversion operations. We also 
emphasize the achievable ’’massive MIMO” effects, even 
for strongly reduced-complexity linear detection techniques, 
provided that the number of BS antennas is much higher 
than the number of antennas which are jointly employed in 
the terminals of the multiple autonomous users. 

1. Introduction 

Cyclic-Prefix (CP)-assisted block transmission schemes 
were proposed and developed, in the last two decades, 
for broadband wireless systems, which have to deal 
with strongly frequency-selective fading channel condi¬ 
tions. These schemes take advantage of current low- 
cost, flexible, Fast Fourier Transform (FFT)-based signal 
processing technology, with both Orthogonal Frequency 
Division Multiplexing (OFDM) and SC/FDE alternative 
choices [1], [2], [3]. Mixed air interface solutions, with 
OFDM for the downlink and SC/FDE for the uplink, 
as proposed in [2], are now widely accepted; the main 
reason for replacing OFDM by SC/FDE, with regard to 
uplink transmission, is the lower envelope fiuctuation of 
the transmitted signals when data symbols are directly 
defined in the time domain, leading to reduced power 
amplification problems at the mobile terminals. 


The development of MIMO technologies has been 
crucial for the ’’success story” of broadband wireless 
communications in the last two decades. Through spatial 
multiplexing schemes, following and extending ideas early 
presented in [4], MIMO systems are currently able to 
provide very high bandwidth efficiencies and a reliable 
radiotransmission at data rates beyond 1 Gigabit/s. Appro¬ 
priate MIMO detection schemes, offering a range of per¬ 
formance/complexity tradeoffs [5] - and also joint iterative 
detection and decoding schemes [6], have been essential 
for the technological improvements in this area. In the 
last decade, MU-MIMO systems have been successfully 
implemented and introduced in several broadband com¬ 
munication standards [7]; in such ’’space division mutiple 
access” systems, the more antennas the BS is equipped 
with, the more users can simultaneously communicate in 
the same time-frequency resource. 

Recently, the adoption of MU-MIMO systems with a 
very large number of antennas in the BS, much larger 
than the number of mobile terminal (MT) antennas in 
its cell, was proposed in [8]. This ’’massive MIMO” 
approach has been shown to be recommendable for 
several reasons [8], [9], [10]: simple linear processing 
for MIMO detection/precoding (uplink/downlink), namely 
when using OFDM for broadband block transmission, 
becomes nearly optimal; both Multi-User Interference 
(MUI)/Multi-Stream Interference (MSI) effects and fast 
fading effects of multipath propagation tend to disappear; 
both power efficiency and bandwidth efficiency become 
substantially increased. 

This paper deals with SC/FDE for the uplink of a 
MU-MIMO system where the BS is constrained to adopt 
low-complexity detection techniques but can be equipped 
with a large number of receiver antennas. In this context, 
either a linear detection or a reduced-complexity iterative 
DF detection are considered. As to the linear detection 
alternative, we include both the optimum Minimum Mean- 
Squared Error (MMSE) [11] and the quite simple Matched 
Filter (ME) detection cases. The iterative DF detection 
alternative, which resorts to joint cancellation of estimated 
MUI/MSI and Inter-Symbol Interference (ISI), does not 
involve channel decoding,differently from the iterative 
receiver technique of [6]; it can be regarded as an exten- 



sion to the multi-input context of the reduced-complexity 
iterative receiver techniques previously considered for 
SIMO systems by the authors (see [12], [13], [14] and 
the references therein). 

Regarding performance evaluation by simulation, ap¬ 
propriate semi-analytical methods are proposed, combin¬ 
ing simulated channel realizations and analytical com¬ 
putations of BER performance which are conditional on 
those channel realizations; selected analytical and semi- 
analytical performance bounds and a simple characteriza¬ 
tion of ’’massive MIMO” effects are also provided. This 
paper shows and discusses a set of numerical performance 
results. The main conclusions of the paper are presented 
in the final section. 


II. System Model 


A. SC/FDE for MU-MIMO Uplink Block Transmission 

We consider here a CP-assisted SC/FDE block transmis¬ 
sion, within a MU-MIMO system with Nt TX antennas 
and Nr RX antennas; for example, but not necessarily, 
one antenna per Mobile Terminal (MT). We assume, 
in the jth TX antenna (j = 1,2,...,A/t) a length-A/ 


block s 


(j) = 




of time-domain data 


symbols in accordance with the corresponding binary 
data block and the selected 4-QAM constellation under 
a Gray mapping rule. The insertion of a length-Lg CP, 
long enough to cope with the time-dispersive effects of 
multipath propagation, is also assumed. 

By using the frequency-domain version of 

the time-domain data block given by 
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{j = 1 , 2 ,-•• ,Nt), we can describe the frequency- 
domain transmission rule as follows, for any subchannel 
k {k = 0,lr- ,N-1): 
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channel realization, and = 
is the resulting, frequency-domain, ’’output vector’ 

As to a given MIMO channel realization, it should 
be noted that the Channel Frequency Response (CFR) 

Udd) = , concerning the an¬ 

tenna pair {i,j), is the DFT of the Channel Impulse 
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= 0 for n > Ls {n = 0,1,..., A/ — 1). Regarding a 
statistical channel model - which encompasses all possible 
channel realizations -, let us assume that E = 0 


and E 
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(of course, with = 0 for n > Lg), the average bit 

energy at each BS antenna is given by 
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B. Linear Detection Techniques 

An appropriate linear detector can be implemented 
by resorting to frequency-domain processing. After CP 
removal, a DFT operation leads to the required set 
{Y,; fc = 0 ,l,... ,N — 1} of length-A/i^ inputs to the 
frequency-domain detector (Yk given by Q); it works, 
for each k, as shown in Fig. leading to a set 

|Y/c; /c = 0, !,••• ,A/ — ijof length-A /t outputs Y^ = 

[Yd\Yd\--- yNf (fc = 0,l,-- - ,N-l). 

When Nt < Nr, possibly with Nr ^ 1, either an 
MMSE, frequency-domain, optimum linear detection or 
a reduced-complexity, frequency-domain, linear detection 
can be considered. In all cases, the detection matrix, for 
each subchannel k (/c = 0, 1 ,..., A/ — 1 ) can be written as 


Dz. 




(4) 


where is the conjugate transpose of the estimated MU- 
MIMO channel matrix and is a selected Nr x Nr 
matrix, possibly depending on H/.. Therefore, Y^ = 
DfeYfc = Yk Sit the output of the frequency- 

domain linear detector (see Fig. ^sl)). 

It should be noted that the jth component of H^Y/c 

is given by Y, Hd^’ Y^ {j = 1,2, - ■■ ,Nt)- this 

means that the H;. factor provides Nr Maximal Ratio 
Combining (MRC) procedures, one per MT antenna, all 
of them based on an appropriate MF for each component 
of the length-A/i^ received vector at subchannel k. 

For a MMSE detection - the optimum linear detection 
- or a Zero Forcing (ZF) detection [5], [11], an inversion 
of each Nr x Nr Ak matrix (/c = 0,1, • • • ,N — 1) is 
required. 

A reduced-complexity linear detection can be achieved 
by using an Nr x Nr diagonal matrix A^. The easiest 
implementation corresponds to adopting an identity matrix 
Ak = Iat^. Of course, Y>k = and Yk = H^Y/. 
when Ak = Iat^, which means an ”MF detection”, 
actually not requiring a matrix inversion. 

For a given channel realization H/. and a given detec¬ 
tion matrix D/^, which depends on the estimated channel 











































(Matched Filtering/Maximal Ratio Combining) 

(a) 



(b) 

Fig. 1. Frequency-domain linear detection (a) and iterative DF detection 
combining a linear MF detection and interference cancellation in the 
frequency domain (b). 


realization H/., the output of the frequency-domain detec¬ 
tor is given by 


Y/c = D/cY/c = Tj^Sk + N'/g, (5) 

where Tk = and N'k = D/eN/e. 

With SC/FDE (time-domain data symbols), an Inverse 
Discrete Fourier Transform (IDFT) is required for each 
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vector. The nth com- 


ponent of the resulting length-Y IDFT 
vector can be written as 


C Low-Complexity Iterative DF Technique 

A low-complexity iterative DF technique can be easily 
devised having in mind eq. 0 . This frequency-domain 
nonlinear detection technique combines the use of a linear 
detector and, for all iterations after the initial iteration 
(i.e., for p > 1), a cancellation of residual MUI - and 
residual MSI, when some users adopt several TX antennas 
for spatial multiplexing purposes - as well as residual ISI; 
such cancellation is based on the estimated data block 
which is provided by the preceding iteration and fed 
back to the frequency-domain detector. The output of this 
frequency-domain detector, for iteration p, is as follows: 


Y'M = Yfe + 


7(p) - ffc(p) 


Sfe(p- 1), 


(9) 


[k = ,N -l;p> 1 (for p= 1, Y^(p) = Y*) , 

where ffe(p) = Dfe(p)Hfc - with Dfe(p) denoting 
the detection matrix employed in iteration p - 
and the entries of the diagonal matrix 


7 (p) are given by y^\p) = ^ 
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The implementation of this iterative DF technique is 
especially simple when T>k{p) = for any p, i.e., when 
a linear MF detector is adopted as shown in Fig. (b) for 
all iterations; the matrix inversion which is inherent to 
more sophisticated linear detectors is then avoided. On 
the other hand, a slightly improved performance can be 
achieved through an increased implementation complexity, 
by feeding back vectors of soft (instead of hard) time- 
domain symbol decisions for interference cancellation. 


III. Evaluation of the Achievable Detection 
Performances 


^i) = jU)sU) + ISI + MUI + 'Gaussian noise', ( 6 ) 
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Therefore, Yk can be written as 


Y/c = 7 S/C + (F/c — 7 ) S/c + D/jN/c, (7) 


where 7 is a diagonal Nt x Nt matrix with (j, j) entries 
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the four terms in the right-hand side of eq. 0 are 
concerned, respectively, to ’’signal”, ISI, MUI/MSI and 
’’Gaussian noise”, at subchannel k. 


A. Semi-analytical Performance Evaluation 

Regarding evaluation of detection performances by 
simulation, simple semi-analytical methods are presented 
here, for the detection techniques of subsecs |II-B| and 
|II-C[ both combining simulated channel realizations and 
analytical computations of BER performance which are 
conditional on those channel realizations. In all cases, the 
conditional BER values are directly computed by resort¬ 
ing to a Signal-to-Interference-plus-Noise-Ratio (SINR), 
under the assumption that the ’’interference” has a quasi- 
Gaussian nature. These ratios are simply derived in accor¬ 
dance with the channel realization (/c = 0 ,l,-' - ,Y — 
1). Of course, for concluding the BER computation in each 
case - involving random generation of a large number of 
channel realizations and conditional BER computations 
- a complementary averaging operation over the set of 
channel realizations is required. 

When using a linear detection technique (sec. II-B[ ), it 
is easy to conclude, having in mind that the ’’signal- 
to-interference-plus-noise” ratio concerning the jth input 
of the MU-MIMO system is given by 










































TABLE I 

SIMO (1 X Nr) Performance Bounds. 
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For 4-QAM transmission, the resulting BERj (j = 
1,2,-•• ,A^t) - conditional on the channel realization 
{H/c; /c = 0,1, • • • , — 1} - is given by 


BERn 




(11) 


(where Q{.) is the Gaussian error function) with SINRj 

Nt 

as computed above, and BER = 7 ^ E BERj. 

^ i-i 

When using the iterative DF technique of Fig. (b), 
eq. replaced by BERj{p) ^ Q (^^ySINRj (p)^ 

for p> 1; certainly, D/c(p) = Hf and Tk{p) = Hf 
in the computation of SINRj{p). By assuming a perfect 
channel estimation , this computation for p > 

1 can be made especially simple, taking into account that 
Y'{p) = 7 Sfe + {Tk - 7 ) (Sfe - Skip - 1)) + DfeNfe 
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[11]. It can be shown 


II-B 


Successively improved performance bounds can be 
obtained as follows (see Table |^: also under D/^ = 
H- qTatt.) ^ H^, with the same Nr but Nr = 1, 
which corresponds to a SIMO/Linear Detection Bound 
(LDB); under T>k = H^, with Nt = 1 and the same Nr, 
by suppressing the resulting first term (ISI) in the denom¬ 
inator of eq. which corresponds to a SIMO/MFB; 
when replacing the fading channel by an AWGN channel 
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for each (i, 1) antenna pair (i = 1,2,-•• ,Nr) in the 
SIMO/MFB, which corresponds to a SIMO/AWGN/MFB. 

It should be noted that the BER curve for the 
SIMO/AWGN/MFB actually corresponds to the achiev¬ 
able BER performance (against White Gaussian Noise) in 
a SIMO 1 X Nr system with single-path propagation for 
all (i, 1) antenna pairs, provided that an ME detection, 
under ideal channel estimation, is adopted: 


BER = qU2„N^^^ 


(13) 


C. Massive MIMO effects 

When Nr > Nr, both the MUI/MSI effects and the 
effects of multipath propagation (fading, ISI) tend to dis¬ 
appear: consequently, the BER performances for the MU- 
MIMO Nt X Nr Rayleigh fading channel become very 
close to those concerning a SIMO 1 x Nr channel with 
single-path propagation for all Nr TX/RX antenna pairs. 
The achievable performances under a ’’truly massive” 
MU-MIMO implementation can be analytically derived 
as explained in the following. 

Entries of are i.i.d. Gaussian-distributed random 
variables with zero mean and variance P^- Therefore, 
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B. Reference MMSE Performance and SIMO Performance 
Bounds 

When adopting an ”MMSE detector” and a 
perfect channel estimation is assumed, D/. = 


lim 
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Consequently, for Nr ^ Nr, 
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that the resulting SINRj - which can be used for 
computing BERj, and then BER - can be wri tten as 
SINRj = defined in subsec. 
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Fig. 2. BER performances for SC/FDE-based MU-MIMO, with Nt = 10 and Nji = 30 (a), 50 (b) or 100 (c), when using linear detection (MF 
and MMSE) [with the SIMO/LDB (dashed line), the SIMO/MFB and the SIMO/AWGN/MFB (lx Nji) performances also included]. 





Fig. 3. BER performances, at the several iterations (p = 1, 2, 3, 4) for SC/FDE-based MU-MIMO, with Nj- = 10 and Nji = 30 (a), 50 (b) or 
100 (c), when using the iterative DF detection technique. 


(by assuming that T>k = H^). When Nr > Nr, 

SINK, =» N„ (18) 

= 2r]NR§^, 

which implies that BER ^ Q ^^2riNR^\ 
i.e. a BER performance closely approximating tne 
SIMO/AWGN/MFB (eq. ([^). 

When Nr ^ Nr, it should also be noted - having in 
mind the equations ini and ( p^ - that the MMSE linear 
detection becomes practically equivalent to a ME linear 
detection, since 

Dfc = (HfHfc + alA(^)“'Hf (19) 

with /3 = (assuming a perfect channel estima¬ 

tion). Of course, the corresponding matrix is then 
r, « leading to 7 ^^'^ « 

~ a+NnK ^ ,Nt). Therefore, the 

resulting SINR’s are as expected, under a perfect chan- 
nel estimation {siNRj = 

Clearly, when Nr > Nr, the MUI/MSI effects as 
well as both the fading and the ISI effects of multipath 
propagation become vanishingly small, leading to a close 
approximation to the SIMO/AWGN/MEB reference per¬ 
formance. 

IV. Numerical Results and Discussion 

The set of performance results which are presented here 
are concerned to SC/EDE uplink block transmission, with 


N = 256 and Ls = 64 in a MU-MIMO Nr x Nr 
Rayleigh fading channel. Perfect channel estimation and 
perfect power control are assumed. The fading effects 
regarding the several TX/RX antenna pairs are supposed 
to be uncorrelated, with independent zero-mean complex 
Gaussian coefficients assumed to have variances 

= 1 — n = 0,1,...,63 (Pn = 0 for n = 
64,65,..., 255). 

The accuracy of performance results obtained through 
the semi-analytical simulation methods of sec. |I^ was 
assessed by means of parallel conventional Monte Carlo 
simulations (involving an error counting procedure), 
which correspond to the superposed dots in the several 
BER performance curves of Eigs.[^and[^ concerning the 
MU-MIMO system. 

Eig. shows the simulated BER performances for an 
SC/EDE-based MU-MIMO uplink and three possibilities 
regarding Nr for Nr = 10, when using two linear de¬ 
tection techniques: optimum (MMSE) detection; reduced- 
complexity (ME) detection. In each subfigure, for the 
sake of comparisons, we also include the SIMO 1 x Nr 
performance bounds of Table |T| In the simulation results 
concerning each subfigure of Eig. the five BER perfor¬ 
mance curves are ordered, from the worst to the best, as 
follows: Nr X Nr MU-MIMO with reduced-complexity 
(ME) linear detection; Nr x Nr MU-MIMO with MMSE 
detection; SIMO/LDB (1 x Nr); SIMO/MEB (1 x Nr); 
SIMO/AWGN/MEB (1 x Nr) [practically superposed to 
the SIMO/MEB curve]. These results clearly show that the 
















performance degradation which is inherent to the reduced- 
complexity linear detection technique (MF) - as compared 
with the MMSE linear detection - can be made quite 
small, by increasing Nr significantly; they also show 
that, under highly increased Nr values, the ’’MUI/MSI- 
free” SIMO (multipath) performance and the ultimate 
bound - the ”MUI/MSI-free and ISI & fading-free” SIMO 
(single-path) performance - can be closely approximated, 
even when adopting the reduced-complexity linear detec¬ 
tion. This figure emphasizes a ’’massive MIMO” effect 
when Nr ^ Nr, which leads to BER performances 
very close to the ultimate ”MUI/MSI-free and ISI & 
fading-free” SIMO (single-path) performance bound (the 
SIMO/AWGN/MFB of Table |I]). 

Fig. shows the simulated BER performances for an 
SC/FDE-based MU-MIMO uplink and three possibilities 
regarding Nr for Nr = 10, when using the reduced- 
complexity iterative DF detection technique of Fig. (b), 
which does not require matrix inversions. By comparing 
these results to the results of Fig. also for Nr = 10, 
we can conclude that, whenever Nr > 3Nr, the reduced- 
complexity DF technique of Fig. (b) is able to provide 
BER performances which are better than those of the 
MMSE linear detector [while avoiding the inversion of 
a 10 X 10 matrix Ak {k = 0,1,...,A^ — 1)], closely 
approximating the (practically identical) SIMO/MFB and 
the SIMO/AWGN/MFB reference performances after a 
small number of iterations. 



Fig. 4. Uplink transmission scenario with six users and Nr = 10 TX 
antennas. 


It should be noted that the numerical results reported 
above, for Nr = 10 TX antennas, are compatible with 
an uplink transmission scenario involving up to 10 users 
(for example, 10 users with one TX antenna per user); a 
specific scenario with Nr = 10 TX antennas, involving 
six users, is depicted in Fig. 

V. Conclusions 

This paper was dedicated to the uplink detection and 
performance evaluation for a MU-MIMO system with 
SC/FDE transmission, when adopting a large number of 
antennas and low-complexity detection techniques at the 
BS. With the help of selected numerical performance 
results, discussed in detail in Section ||Vj we show that 
a moderately large number of BS antennas (say. Nr = 
3Nr) is enough to closely approximate the SIMO/MFB 
performance - and also the SIMO/ AWGN/MF B perfor¬ 
mance, expressed as BER = Q (^^J2r]NR^^ -, espe¬ 
cially when using the suggested low-complexity iterative 
DF technique, which does not require Nr x Nr matrix 
inversion. We also emphasize the ’’massive MIMO” effects 
provided by a number of BS antennas much higher than 
the number of antennas which are jointly employed in the 
terminals of the multiple autonomous users, even when 
strongly reduced-complexity linear detection techniques - 
such as the so-called ”MF detection -, are adopted. 

The accuracy of performance results obtained by semi- 
analytical means, much less time-consuming than conven¬ 
tional, ’error counting’-based, Monte Carlo simulations - 
was also demonstrated. The proposed performance eval¬ 
uation method can be very useful for rapidly knowing 
”how many antennas do we need in the BS?”, for a 
given number of antennas jointly employed in the user 
terminals. 
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